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ABSTRACT 



Aims. The binary pulsar PSRJ181 1-1736 has been identified, since its discovery, as a member of a double neutron star system. 
Observations of such binary pulsars allow the measurement of general relativistic effects, which in turn lead to information about the 
orbiting objects, to their binary evolution and, in a few cases, to tests of theories of gravity. 

Methods. Regular timing observations have since 2000 been carried out with three of the largest European radio telescopes involved 
in pulsar research. Pulse times of arrival were determined by convolving the observed profiles with standard templates, and were fitted 
to a model that takes into account general relativistic effects in binary systems. The prospects of continued observations were studied 
with simulated timing data. Pulse scattering times were measured using dedicated observations at 1.4 GHz and at 3.1 GHz, and the 
corresponding spectral index has also been determined. The possibility of detecting the yet unseen companion as a radio pulsar was 
investigated as a function of pulse period, observing frequency and flux density of the source. A study of the natal kick received by 
the younger neutron star at birth was performed considering the total energy and total angular momentum for a two body system. 
Results. We present an up to date and improved timing solution for the binary pulsar PSRJ181 1-1736. One post-Keplerian parameter, 
the relativistic periastron advance, is measured and leads to the determination of the total mass of this binary system. Measured and 
derived parameters strongly support the double neutron star scenario for this system. The pulse profile at 1.4 GHz is heavily broadened 
by interstellar scattering, limiting the timing precision achievable at this frequency. We show that a better precision can be obtained 
with observations at higher frequencies. This would allow one to measure a second Post-Keplerian parameter within a few years. We 
find that interstellar scattering is unlikely to be the reason for the continued failure to detect radio pulsations from the companion of 
PSRJ181 1-1736. The probability distribution that we derive for the amplitude of the kick imparted on the companion neutron star at 
its birth indicates that the kick has been of low amplitude. 

Key words. Pulsar: General, PSRJ181 1-1736 



1. Introduction 

The pulsar J181 1-1736, discover ed during observations of the 
Parkes Multibeam Pulsar Survey) Manchest er et all 12001 1) , has a 
spin period of 104 ms and is member of a 18.8-d, highly eccen- 
tric binary s ystem (|Lvne et all 120001) with an as yet undetected 
companion (lMignanill2000l) . The characteristic age and the esti- 
mated surface magnetic field strength indicate that the pulsar is 
mildly recycled. In such a system, it is expected that the observed 
pulsar is born first in a supernovae (SN) explosion before it un- 
dergoes mass accretion from a high-mass non-degenerate binary 
companion. Parameters which have been measured and derived 
from the best-fit timing solution indicate that the companion is 
quite massive. All thes e elements suggest that the companion is 
also a neutron star (e.g. Bhattacharya & van den Heuvelf T991). 

The conclusion that PSRJ181 1-1736 is a member of the 
small sample of kn own double neutr on star (DNS) systems was 
already reached bv lLvne et al.l d2000t) . Their conclusion was sup- 
ported by a constraint on the total system mass, assuming that 
the observed advance of periastron was totally due to relativistic 
effects. However, the rather long period of the pulsar, combined 
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with the effects of interstellar scattering, resulting in significant 
broadening of the pul se profile a t 1 .4 G Hz, as well as the short 
data span available to lLvne et al.l d2000l) . limited the timing pre- 
cision and hence the accuracy of the total mass measurement. 

Among all DNS systems, PSRJ181 1-1736 has by far the 
longest spin period, the longest orbital period and the highest ec- 
centricity. This may suggest that the evolution of this DNS has 
been different, at least in part, from all other DNS systems. On 
the other hand, PSRJ181 1-1736 well fits the spin period versus 
eccentricity relation for DNS systems (McLaughlin et al. 2005, 
Faulkner et al. 2005). This relation can be simply explained in 
terms of the different lengths of time the pulsar underwent ac- 
cretion, which in turn is related to the mass of the companion 
star before the SN explosion. Moreover, numerical simulations 
dDewi et al 1 120051) show that the spin period versus eccentricity 
relation is recovered assuming that the second born neutron star 
received a low velocity kick at its birth. 

In this paper we report on new timing observations which 
significantly improve on the previously published results. We 
present a study of the observed interstellar scattering and con- 
sider its consequences on the detectability of radio pulses from 
the companion. Finally, we investigate the likely kick velocity 
imparted to the second-born neutron star during its birth in the 
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system's second SN. Observations were carried out as part of a 
coordinated effort using three of the largest steerable radio tele- 
scopes in the world for pulsar timing observations, i.e. the 100- 
m radio telescope at Effelsberg, the 94-m equivalent Westerbork 
Synthesis Radio Telescope (WSRT) and the 76-m Lovell tele- 
scope at Jodrell Bank telescope. This paper is the first in a se- 
ries detailing results of these efforts in establishing a European 
Pulsar Timing Array (EPTA). 

2. Observations 

The binary pulsar J181 1-1736 is one of the binary pulsars reg- 
ularly observed by the EPTA. The aims and objectives of this 
European collaboration include the detection of a cosmological 
gravitational wave background and the project will be described 
in detail in a forthcoming publication. Here we summarize the 
observing systems used while further details can be found in the 
references below. 

2.1. Effelsberg timing 

We made regular timing observations of PSRJ181 1-1736 since 
October 1999 using the 100-m radio telescope of the Max Planck 
Institut filr Radioastronomie in Effelsberg near Bonn. The typ- 
ical observing rate was of 1 observation every two months. An 
overall root-mean-square (RMS) of 538 (is is achieved after ap- 
plying the final timing model. The data were obtained with a 
1.3-1.7 GHz tunable HEMT receiver installed in the primary 
focus of the telescope. The noise temperature of this system is 
25 K, resulting in a system temperature from 30 to 50 K on cold 
sky depending on elevation. The antenna gain at these frequen- 
cies is 1.5 K Jy _1 . 

An intermediate frequency (IF) centred on 150 MHz for 
left-hand (LHC) and right-hand (RHC) circularly polarised sig- 
nals was obtained after down-conversion from a central RF fre- 
quency of usually 1410 MHz. The signals received from the tele- 
scope were acquired and processed with the Effelsberg-Berkeley 
Pulsar Processor (EBPP) which removes the dispersive effects of 
the interstellar med i um o n-line using "coherent de-dispersion" 
dHankins & Rickettl 119751) . Before entering the EBPP, the two 
LHC and RHC signals of the IF are converted to an internal IF 
of 440 MHz. A maximum bandwidth of 2 x 32 x 0.7 MHz = 2 x 
32 MHz was available for the chosen observing frequency and 
DM of the pulsar. It was split into four portions for each of the 
two circular polarisations, which were mixed down to baseband. 
Each portion was then sub-d ivided into eight narrow channels 
via a set of digital filters (Backe r~et all 1 1997b . The outputs of 
each channel were fed into de-disperser boards for coherent on- 
line de-dispersion. In total 64 output signals were detected and 
integrated in phase with the predicted topocentric pulse period. 

A pulse time-of-arrival (TOA) was calculated for each av- 
erage profile obtained during a 5-10 min observation. During 
this process, the observed time-stamped profile was compared 
to a synthetic template, which was constructed out of 5 Gaussian 
components fit ted to a hig h signal-to-noise standard profile (see 
Kramer et al., 1998I I1999I) . This template matching was d one by 
a leas t-squares fitting of the Fourier-transformed data dTavlorl 
|1992|) . Using the measured time delay between the actual profile 
and the template, the accurate time stamp of the data provided 
by a local H-MASER and corrected off-line to UTC(NIST) 
using recorded information from the satellites of the Global 
Positioning System (GPS), the final TOA was obtained. The un- 
ce rtainty of each TOA was e stim ated using a m ethod described 
by lDowns & Reichleyl d 19831) and|LangeJ d 19991) . 



2.2. Jodrell Bank timing 

Observations of PSRJ181 1-1736 were made regularly using the 
76-m Lovell telesc ope at Jodrell Bank since its discovery in 1997 
dLyne et al.Ll2000h . The typical observing rate was of about 2 ob- 
servations each week, with an overall RMS of 1300/zs after ap- 
plying the final timing model. A cryogenic receiver at 1404 MHz 
was used, and both LHC and RHC signals were observed using 
a 2 x 32 x 1.0-MHz filter bank at 1404 MHz. After detection, 
the signals from the two polarizations were filtered, digitised 
at appropriate sampling intervals, incoherently de-dispersed in 
hardware before being folded on-line with the topocentric pulse 
period and written to disk. Each integration was typically of 1-3 
minutes duration; 6 or 12 such integrations constituted a typi- 
cal observation. Off-line, the profiles were added in polarisation 
pairs before being summed to produce a single total-intensity 
profile. A standard pulse template was fitted to the observed pro- 
files at each frequency to determine the pulse times-of-arrival 
(TOAs). Details of the observing sy stem and the data r eduction 
scheme can be found elsewhere (e.g. lHobbs et alj [2004). 

2.3. Westerbork timing 

Observations of PSRJ181 1-1736 were carried out approxi- 
mately monthly since 1999 August 1st, obtaining an overall tim- 
ing RMS of 659 /vs after applying the final timing model, at a 
central frequency of 1380 MHz and a bandwidth of 80 MHz. 
The two linear polarisations from all 14 telescopes were added 
together in phase by taking account of the relative geometrical 
and instrumental phase dela ys between t hem a nd then passed 
to the PuMa pulsar backend dVoute et all l2002h . The data were 
obtained with the L-band receiver installed in the primary fo- 
cus of the telescopes. The noise temperature of this system is 
25 K, resulting in a system temperature from 30 to 50 K on 
cold sky depending on elevation. The antenna gain at these fre- 
quencies is 1.2 K Jy _1 . PuMa was used in its digital filterbank 
mode whereby the Nyquist sampled signals are Fourier trans- 
formed and the polarisations combined to produce total inten- 
sity (Stokes I) spectra with a total of 512 channels. These spec- 
tra were summed online to give a final sampling time of 409.6 
jus and recorded to hard disk. These spectra were subsequently 
dedispersed and folded with the topocentric period off-line to 
form integrations of a few minutes. TOAs were calculated for 
each profile following a scheme similar to that outlined above 
for Effelsberg data, except a high signal-to-noise standard profile 
was used instead of Gaussian components. In the future, EPTA 
timing will employ an identical synthetic template for all tele- 
scopes. 



3. Data analysis 

The TOAs, corrected to UTC(NIST) via GPS and weighted 
by their individual uncertainties determined in the fitting 
process, were analysed with the TEMPO software package 
dTavlor & Weisbergl 1 19891) . using the DE405 ephemeris of the 
Jet Propulsion Laboratory dStandishl H990I) . TEMPO minimizes 
the sum of weighted squared timing residuals, i.e. the difference 
between observed and model TOAs, yielding a set of improved 
pulsar parameters and post-fit timing residuals. A summary of 
the basic characteristics of each dataset is shown in Table 1 . 

Although the templates used for the three telescope data dif- 
fered, resulting offsets were absorbed in a global least-squares 
fit. Remaining uncertainties were smaller than the typical mea- 
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Table 1. Data sets' characteristics 





Jodrell Bank 


Effelsberg 


Westerbork 


N. of ToAs 


348 


74 


213 


Time Span 


50842-53624 


51490-53624 


51391-53546 


R.M.S. 


1300 


538 


659 




o looo gooo :jooo 

Da Lc (MJD since 50800) 



Fig. 1. Timing residuals after jointly applying the final model to 
all three data sets. Vertical bars represent the ToA's uncertainty. 

surement accuracy of the Jodrell Bank timing data of about 9 

Before all TOAs were combined, preliminary fits were per- 
formed on each dataset alone, in order to study possible system- 
atic differences between the datasets. We applied a small quadra- 
ture addition and a scaling factor to the uncertainties to obtain the 
expected value of a reduced x 2 — 1 f° r eacn dataset. The final 
joint fit to all TOAs resulted in a x 2 value of unity, avoiding the 
need to add further systematic uncertainties. 

Table 2 summarizes all observed timing and some derived 
parameters. For the observed parameters the quoted errors are 
twice the nominal TEMPO errors. For the derived parameters, 
the given uncertainties are computed accordingly. 

The joint fit allowed us to determine the spin, positional and 
Keplerian orbital parameters plus one post-Keplerian parameter 
with a precision better than the best determination from a single 
data set alone. However, the high degree of interstellar scattering 
(Fig. 3) means that further post-Keplerian parameters will be dif- 
ficult to measure with continued observations at this frequency. 
We will discuss the future prospects for higher frequency obser- 
vations in § 5. 

4. The nature of the companion 

In their discovery paper, iLvne et all. J2000) proposed that this 
syste m is a member o f the small class of DNS binaries. Soon 
after, Mignani (2000) reported on optical observations of the 
region surrounding the pulsar position to search for emission 
from the pulsar companion. They detected no emission coinci- 



Table 2. Timing and derived parameters 



Timing parameters 


Joint data sets 


RA (J2000, hh:mm:ss) 


18:11:55.034(3) 


DECL (J2000, deg:mm:ss) 


-17:36:37.7(4) 


Period, P (s) 


0.1041819547968(4) 


Period derivative, P (10~ 19 s s~') 


9.01(5) 


Dispersion Measure, DM (pc cm ^) 


476(5) 


Projected semi-major axis", a sin / (s) 


34.7827(5) 


Eccentricity, e 


0.828011(9) 


Epoch of periastron, To (MJD) 


50875.02452(3) 


Orbital period, Pg (d) 


1 S 770 1 fid 1 {A\ 


Longitude of periastron, ll> (deg) 


127.6577(11) 


Advance of periastron, io (deg yr -1 ) 


0.0090(2) 


I — '1 j ■ . j ~t 1 f\f\ TV JI I ( ' / T \ 

Flux density at 3 100 MHz, S 3100 (mJy) 


0.34(7) 


lime Span (MJD) 


5 1)542-3 J 024 


JN. Ot loAS 




Klvlo \pLS ) 


□ ci 1-70 


Derived parameters* 




Characteristic Age, t c (10 9 yr) 


1.83 


Surface magnetic field, B (10 9 G) 


9.80 


Total Mass M T0T (M Q ) 


2.57(10) 


Mass Function f(M c ) (M Q ) 


0.128121(5) 


Orbital separation A (Is) 


94.4(6) 


Minimum companion mass Mc, m i n (M Q ) 


0.93 



" The projected semi-major axis a sin i is the semi-major axis of the 
projection of the orbit of the pulsar, around the system's center of 
mass, onto the plane containing the line of sight and the line of 
nodes. 

b Characteristic age and surface magnetic field have been calculated 
using standard formulas, namely r c = P/2P and B = 3.2 x 

10 19 Vpp G. The total mass Mtot has been calculated from the rela- 
tivistic periastron advance and the measured Keplerian parameters, 
assuming the validity of general relativity. The minimum compan- 
ion mass was estimated using the observed mass function f(M c ) 
and the lower limit for the total mass, as given by its uncertainty, in 
the case of sin i = 1. For details see Lorimer & Kramer (2005). 



dent with the pulsar position and while not conclusive, the lack 
of emission is at least consistent with the neutron star hypothesis 
for the nature of the companion. 

The derived values for the characteristic age (r c = 1.83 x 10 9 
yrs) and the surface magnetic field (B = 9.8 x 10 9 G), as well 
as the combined values of the spin period (P = 104 ms) and its 
period derivative (P = 9 x IO 19 ), indicate that PSRJ181 1-1736 
is a neutron star that experienced a spin-up phase via accre- 
tion from mass overflowing from its companion. These pa- 
rameters, in conjunction with the measured orbital eccentric- 
ity (e = 0.828), indeed suggest that PSRJ181 1-1736 is a 
mildly recycled pulsar whose companion star was massive 
enough to also undergo a SN explosion. This second SN im- 
parted the actually observed large ecce ntricity to the system (e.g. 
Bhattacharya & van den HeuvellfT99l . 

Our new measurement of the relativistic periastron advance, 
io = 0.0090 + 0.0002 deg yr -1 , allows us to determine the 
value of 2.57+0.10 M for the total mass of the system, assum- 
ing that general relativity is the correct theory of gravity and 
that the observed value is fully due to relativistic effects (e.g. 
iDamour & Deruellall986l) . This value, combined with the mea- 
sured mass function, implies a minimum companion mass of 
0.93 M . 
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Fig. 2. Mass-mass diagram for the binary system hosting 
PSRJ181 1-1736. The shaded area below the dotted curved line 
is excluded because of the geometrical constraint sin i < 1, while 
the area outside the diagonal stripe is excluded by the measure- 
ment of the relativistic periastron advance and the derived value 
for the total mass for this system. 



The value for the total mass is relatively lo w, but very similar 
to the total mass of the double pulsar system dLvne et all [2004) 
and the recently d iscovered DNS system PSR J1756-2251 
dFaulkner et all 120051) . In fact, these systems have neutron star 
companions that have the lowest neutron star masses observed 
so far, M c = 1.25M and M c = 1.1 7M Q (for a recent re- 
view see Stairs 2004), respectively. In Figure [2] we show the 
so-called mass-mass diagram where the pulsar and companion 
masses can be directly compared. The measured value for the 
advance of periastron means that the sum of the masses must 
lie along the diagonal line, while the constraint on the incli- 
nation sin i < 1 excludes the hatched region below the dotted 
line. Assuming that the neutron stars in this system must have 
a mass which is larger than the lowest mass so-far measured, 
i.e. 1.17M , we find that they both have masses in the interval 
1.17 M Q < M P ,M C < 1.5OM . This interval contains all but 
the heaviest neutron stars masses for which a reliable determina- 
tion has been obtained. Using this mass constraint, we can also 
translate this range into lower and upper limits on the inclination 
of the system, i.e. 44 deg i i i50 deg. 

Alternatively, if either the pulsar or the companion have a 
mass equal to the observed median neutron star mass of 1.35 
M dStairsl |2004|) . the other neutron star would have a mass of 
1.22 + O.IOMq. This value is consistent with the lower limit in 
the previous discussion, but this also allows for the possibility 
that one of the two neutron stars has a mass as low as 1.12 M . 

5. Future potential of timing observations 

In order to determine the companion mass without ambiguities, 
it is necessar y to measure a second post-Keplerian (PK) pa- 
rameter (e.g. Damour & Taylorll 1992b . We have investigated the 
possibility of measuring the PK parameter y, which describes 



the combined effect of gravitational redshift and a second or- 
der Doppler effect. For a companion mass of 1.35M , the ex- 
pected value is y — 0.021 ms. Using simulated data sets for the 
presently available timing precision, we estimate that a 3<x detec- 
tion for y is achievable after about 4 more years of observation. 
However, in order to obtain a 10% accuracy in mass measure- 
ment by determining y to a similar precision, several decades of 
observations may be needed. 

From similar simulations, we estimate that PK parameters 
like the rate of orbital decay, Pb, or the Shapiro delay, are un- 
measurable in this system, unless significant improvements in 
timing precision can be obtained. For a companion mass of 1 .35 
M Q , Pb is only -9.4 x 10~ 15 s s _1 and we expect an amplitude 
of only 6 fis for the amplitude of the Shapiro dela y. We also note 
that the effects of geodetic precession (see e.g. [Kramer] |l998) 
will not be measurable within a reasonable time, as it has a pe- 
riod of order of 10 5 years. 

6. Improving timing precision 

It is obvious that the measurement of further PK parameters will 
only be possible if higher timing precision can be achieved for 
this pulsar. For instance, if a precision of 50^s could be obtained, 
y could be measurable to a 10% accuracy after a total of just 5 
yr of observations, while a 3cr detection of the orbital decay may 
be achieved after about 6 yr. 

One way to achieve higher timing precision is to detect nar- 
row features in the observed pulse profile by means of higher ef- 
fective time resolution. This is most commonly achieved through 
better correction for dispersion smearing that is caused by the 
radio signal's passage through the ionized interstellar medium. 
While this effect is actually completel y removed by the use co- 
herent de-dispersion techniques (see lHankins & Ricketll 1 19751) 
at some of our telescopes, it is apparent that the current timing 
precision is limited instead b y broadening of th e pulse profile 
due to interstellar scattering dLohmer et al.l l2004. and references 
therein). Indeed, the pulse profile at 1 .4 GHz shows a strong scat- 
tering tail (see Fig. 3) which prevents a highly accurate determi- 
nation of the pulse time of arrival. As scattering is a strong func- 
tion of observing frequency, we can expect to reduce its effect, 
and hence to enable higher timing precision, by using timing ob- 
servations at frequencies above 1 .4 GHz. 

We obtained observations at 3.1 GHz that confirm this ex- 
pectation. The pulse profile obtained at this frequency shows 
no evidence of interstellar scattering, and its width at 10% is 
only 7.3 ms. This is a great improvement with respect to the 
1.4 GHz profile, whose 10% width is 58.3 ms. A flux density of 
0.34+0.07mJy measured at 3. 1 GHz suggests that regular timing 
observations at this frequency should be possible and should sig- 
nificantly improve the achievable timing precision. This would 
allow us to measure a second PK parameter to an accuracy that 
is sufficiently precise to determine the companion mass. 

Using the data available at 1.4 GHz and 3.1 GHz, we com- 
puted spectral indexes for flux density and scattering time. For 
the flux density we obtain a spectral index of yS = —1.8 ± 0.6. 
Subdividing our observing band at 1 .4 GHz we obtain two dif- 
ferent profiles that we use to measure the p ulse scatter timescale 
t by applying the technique described in Loh mer et all 1200 ll 
We convolve the 3.1 GHz-profile, assumed to represent the true 
pulse shape, with an exponential scattering tail and obtain scat- 
tering times by a least-square comparison of the convolved 
profile with the observed pulse shape. At 1.284 GHz we find 
t.5 = 16.9 ms, and t s = 10.6 ms at 1.464 GHz, respec- 
tively. This results in a spectral index a of the scattering time, 



A.Corongiu et al.: The binary pulsar PSRJ181 1-1736. 



5 




0.4 0.6 
Pulse phase 



Fig.3. Pulses' profiles of PSR J1811-1736 at 1.4GHz (bottom 
panel) and 3.1 GHz (top panel). Both profiles have been obtained 
with 10 minutes observations performed with the Parkes radio 
telescope in February 2005. 



i.e. r oc v~°, of a — 3.5 ± 0.1. Such a measured value agrees 
very well with analogous results from Lohmer et al. (120011: 12004) 
who determined a for a number of pulsars with very high disper- 
sion measures. 

The measured spectral index of the scattering time is also 
consistent with the fact that the pulsar has not been detected at 
frequencies below 1 GHz. For example at 400 MHz we calculate 
t s ~ 1 s which is almost an order of magnitude greater than the 
spin period of the pulsar thus making it impossible to detect it as 
a pulsating source. 



7. Previous searches for pulsations from the 
companion 

Searches for pulsations from the binary companion of 
PSR J181 1-1736 have been performed on Parkes and Effelsberg 
data. Parkes observ ations have be e n inv estigated with the proce- 
dure described in Faul kner et alj d2004l) . while Effelsber g data 
have been processed using the procedure described in Klein 
(2004). Both searches were unsuccessful in detecting any evi- 
dence of pulsation. The very high value of the dispersion mea- 
sure for this system may suggest that the interstellar scatter- 
ing is responsible for the failure in detecting any pulsation. 
Therefore we studied the possible impact of this phenomenon on 
our searches for pulsations from the companion of PSRJ1811- 
1736. 

We considered 1 hr observations done with the Effelsberg 
telescope using either the 20 cm (1.4 GHz) or the 11cm 
(2.7 GHz) receiver, exploring a range of possible flux densities 
(S = 0.05,0.5, lmJy) and a detection in a signal-to-noise ra- 
tio threshold of S /N = 10. Using the DM of the observed pulsar, 
and assuming Effelsberg observations at 1 .4 GHz, we find a min- 
imum detectable period of P m j n = 750 ms for a flux density of 
5=50 fiiy, while even for the flux densities of S — 500 /uJy and 



5=1 mJy periods below ~10ms become undetectable at the 
observing frequency of 1 .4 GHz. 

At 2.7 GHz, the system performance of the Effelsberg tele- 
scope allows for an antenna gain of G =1.5 K3y l with a sys- 
tem temperature r sys = 17 K. Using these parameters, we ob- 
tain minimum periods of f m i n = 110 ms, P m [ n = 2.5 ms and 
Anin = 1.6 ms for flux densities S 1 2.7 GHz = 50 yi/Jy, 500 piy and 
1 mJy respectively. 

When searching for pulsations from the binary companion 
of a Galactic recycled pulsar in a double neutron star system, it 
is more likely that the companion is a young pulsar with rather 
ordinary spin parame ters, as found for PSR J0737-3039B in the 
double pulsar system dLvne et al.L 12004). Our lower limits on the 
minimum detectable period therefore suggests that interstellar 
scattering should not have prevented the detection of the com- 
panion, unless it were a very fast spinning or very faint source. 

8. Constraints on the kick velocity of the second SN 
explosion 

The large eccentricity of J181 1-1736 system can be ascribed to 
a sudden loss of mass which results in a change of the orbital pa- 
rameters. Such a sudden loss of mass can be attributed to the SN 
explosion that formed the younger unseen neutron st ar compan- 
ion (see, e.g.,[Bhattacharva & van den Heuvelll 19911) . Under the 
hypothesis of a symmetric explosion, simple calculations show 
that the binary survives this event only if the expelled mass M exp 
is less than half of the total mass of the binary before the ex- 
plosion (pre-SN binary). The induced eccentricity is a simple 
function of the amount of the expelled mass: e - M exp /MjoT, 
where Mtot is the total mass of the pre-SN binary. In the case of 
the binary system hosting PSR J181 1-1736, the measured eccen- 
tricity, e = 0.828, and the derived total mass, M\, m = 2.57 M , 
would imply a total mass Mtot = 4.7 M Q for the pre-SN binary. 

The high space velocities measured for isolated pulsars indi- 
cate that neutron stars may receive a kick when formed, with 
an unpredictable amplitud e and d i rection (Hansen & Phi nnevi 



11997b ICordes & Chernoffl Il998t lArzoumanian et all l2002h . 



Such kicks imparted to the newly formed neutron stars are 
caused by asymmetric supernova explosions. If an asymmetric 
SN explosion occurs in a binary system, the survival and the 
eventual post-SN binary parameters are jointly determined by 
the mass loss and the vector representing the velocity imparted 
to the neutron star. In this case a simple survival condition like 
the one derived for the symmetric explosion case cannot be de- 
termined. 

A correlation between the pulsar's spin period and or- 
bital eccentricity has recently been found for DNS systems 
dMcLaughlin et al.l 120051 iFaulkner et alj 120051) . A numerical 
simulation by iDewi et ail (120051) linked this correlation to the 
typical amplitude o f the kick velocity received by the younger 
neutron star at birth. [Dewi et al.l (12005 ) found that the spin period 
versus eccentricity correlation is recovered if the typical kick 
amplitude satisfies the condition Vk ~ 50kms _1 . 

To investigate the nature of the kick received by the younger 
neutron star in this system we considered as the pre-SN binary 
a binary system containing a neutron star and a helium star. 
We then constrained the total mass of this system by combin- 
ing our results on the total mass o f the actual binary with the 
range given by IDewi &~P ols (2003) for the mass of the helium 
star that was the companion of PSR J1811-1736 before the ex- 
plosion. The helium st a r mas s range 2.8M Q < Mq < 5.OM 
given by IDewi & Polsl d2003l) leads to a total mass range of 
4.0M Q < M TO T < 6.5M Q . 
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Binary parameters for the pre-SN binary have been chosen 
as follows. The eccentricity has been assumed negligible, since 
the accretion phase responsible for spinning up the pulsar also 
provided strong tidal forces that circularised the orbit. The or- 
bital separation has been constrained to be between the the min- 
imum (pericentric) and maximum (apocentric) distance between 
the two neutron stars in the post-SN binary. This statement can 
be justified as follows. The typical velocity of the expelled mat- 
ter in a SN explosion is close to the speed of light, while the typ- 
ical orbital velocity of the stars in a binary system is of order of 
~ 100 km s _I , using for the total mass any value in the range we 
used for Mtot and a value for the orbital separation comparable 
to the one for the present binary system, i.e. few light-minutes 
(see the discussion in the next paragraph of the post-SN binary 
evolution due to general relativistic effects). This means that the 
change in position of the two stars is negligible if compared to 
the change in position of the expelled matter. The time required 
for the binary system to do the transition from the pre-SN to 
the post-SN binary is the time required by the expelled matter 
to travel along a path as long as the orbital separation, i.e. few 
minutes. After such an elapsed time the matter expelled in the 
SN explosion encloses both stars and has no more gravitational 
effects on their binary motion. This time is also much shorter 
than the orbital period of few days for a pre-SN binary like the 
one we are considering. This means that during this transition 
the positions of the two stars remained unchanged, and their dis- 
tance was a distance periodically assumed by the two stars also 
in their orbital motion in the post-SN binary. 

To make a fully consistent comparison between the actually 
observed binary and the eccentric binary that emerged from the 
last SN explosion (post-SN binary) one has to take into account 
the secular changes of the orbital parameters caused by general 
relativistic effects. In order to do this one needs to have an es- 
timate for the time since the last SN. The only timescale that is 
available to us is the characteristic age of the first born pulsar. We 
then find binary parameters that are consistent with the present 
system values within their uncertainties. Given the well known 
uncertainties in this age estimation, we considered the possibil- 
ity that the present binary is in fact up to ten times older than 
suggested by the characteristic age (i.e 1.8xl0 10 yr). Even when 
considering this extreme age, we find that our results remain un- 
affected. We consequently decided to use as post-SN orbital pa- 
rameters the same values we measure today. 

By insisting that the total energy and total angular momen- 
tum, calculated in the center-of-mass frame, of the post-SN bi- 
nary and the present systems are conserved, we can combine 
these terms to obtain an equation for the kick amplitude, as a 
function of the two angles, representing its direction in a suitable 
reference frame, and the total mass and orbital separation before 
the explosion. We assumed that the probability of occurrence for 
any given kick vector is proportional to the solid angle described 
by the direction of the kick in spherical coordinates and then cal- 
culated the probability of having a kick velocity lower than some 
fixed values. We chose the values of 50, 100 and 150 km s _1 ( 
hereafter P50, P100 and P150 respectively.). Figure|4] shows that 
P50 is not negligible if the total mass of the pre-SN binary sys- 
tem is lower than 6M . Moreover, all considered probabilities 
peak in correspondence of a total pre-SN mass of 4.70 M Q , cor- 
responding to the null kick case. These results lead to the con- 
clusion that the younger neutron star in this system received a 
low velocity kick and is thus similar to all other known DNS 
systems, which all have tighter orbits. 

Nevertheless, the binary system containing PSRJ181 1-1736 
is much wider than all other known DNS systems. This may in- 
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Fig. 4. Probabilities to have a kick velocity lower than or equal to 
50 (lower line), 100 and 150 (upper line) km s _1 as a function of 
the pre-SN total mass. All probabilities peak in correspondence 
of a pre-SN total mass Mtot = 4.7OM , which is the mass of the 
binary system before the explosion in the case of a symmetric 
SN. The probability to have a kick velocity lower than 50 km 
s is not negligible for all but the highest considered values for 
the binary pre-SN mass. 



dicate that the binary evolution of this system may have been 
(at least partially) different. In particular the wide orbital sepa- 
ration for this system may be compatible with an evolution dur- 
ing which the pu lsar's progeni t or a voided completely a common 
envelope phase (Dewi & Pols. 12003b or that this phase was too 
short to sufficiently reduce the orbital separation. Moreover if 
the spin-up occurred via the stellar wind of the giant companion 
then the system would ten d to be wider due to the isotropic mass 
loss from the companion dDewi et al 1 12003) . 

9. Summary & Conclusions 

We have presented an improved timing solution for the binary 
pulsar J181 1-1736. This solution improves the previously mea- 
sured values for the spin and Keplerian orbital parameters and 
one post-Keplerian orbital parameter, the periastron advance. 
These results would not have been achieved without data from 
the three telescopes used and are the first obtained as part of the 
European Pulsar Timing Array (EPTA) collaboration. 

The measured values for the spin period and its first deriva- 
tive are typical of a mildly recycled neutron star, while the high 
eccentricity of the binary system can be seen as a signature of the 
SN explosion that interrupted the mass transfer from the com- 
panion to the accreting neutron star. This is likely to have oc- 
curred before the pulsar could reach spin periods typical of the 
fully recycled (i.e. millisecond) pulsars. This leads to the conclu- 
sion that PSRJ181 1-1736 is a member of a DNS binary system. 

The determined value of the periastron advance provides fur- 
ther confirmation for this scenario as it suggests a total mass 
of the system of M, , - 2.57±0.10 M . This value is similar to 
the total mass of two other DNS systems, i.e. the double pulsar 
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dLvne etalll200l and PSR J1756-2251 dFaulkner et alll2005h . 
In both these systems, the non-recycled neutron stars is very 
light. Assuming that PSR J181 1-1736 is a neutron star with a 
mass within the currently measured mass range for neutron stars, 
we find the companion mass to lie in the same range. Using these 
arguments we determine the inclination of the orbital plane to 
within ~ 6 degrees. 

We also investigated the possibility of measuring a second 
post-Keplerian parameter, in order to determine both masses 
and thus to definitively determine the nature of the companion. 
Unfortunately, the pulse profile at 1 .4 GHz is heavily broadened 
by interstellar scattering which limits the timing precision and 
means that a second post-Keplerian parameter is not measur- 
able within a reasonable amount of time with observations at 
that frequency. However we find that at 3 GHz the scattering 
is sufficiently reduced and the flux density is sufficiently high 
that higher precision timing will be possible at this frequency. 
Comparing the pulse profiles at 1 .4 and 3 GHz we find that the 
scattering timescale for this pulsar scales with frequency with a 
power law of index a = 3.5 + 0.1 which is in excellent agree- 
ment with earlier results on high dispersion measure pulsars. 

Considering the effects of the interstellar scattering on the 
detectability of pulsations from the companion, we find that the 
minimum detectable period is longer at lower frequencies and 
for fainter objects. In general, we do not expect interstellar scat- 
tering to be the cause for the continued non-detection of the com- 
panion neutron star. 

The orbital separation for this system is much wider than for 
all other DNS systems, and it suggests that its binary evolution 
has been different. One explanation invokes the lack of a com- 
mon envelope phase, during which the size of the orbit shrinks 
due to tidal f orces in the envelo pe of the companion star. Another 
explanation dDewi et all 120051) . not necessarily conflicting with 
the previous one, invokes a different mass transfer mechanism in 
the spin-up phase of the pulsar, namely via stellar wind, while 
all other recycled pulsars in the known DNS have been spun up 
via Roche lobe overflow mass transfer. 

Finally, we investigated the kick imparted to the second born 
neutron star during the second SN. We find that for realistic val- 
ues of the total mass of the pre-SN binary, the kick velocity 
has a not negligible probability of being lower than 50 km s _1 . 
This constraint is common to all DNS systems, as shown by 
Dew i et all d2005l) . This evidence for a low amplitude asymmet- 
ric kick received by the younger neutron star may be the conse- 
quence of the effects of binary evolution on a star that undergoes 
a SN explosion, effects that somehow are able to tune the ampli- 
tude of such kick. 
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